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SUMMARY 


/ 

This  memorandum  outlines  the  mathematical  formulation  of  a 
polarimetric  theory  for  radar  scattering.  The  emphasis  is  placed 
on  physical  interpretation  of  some  fundamental  results  from  the 
theory  of  nonsingular  linear  transformations  and  the  general 
scattering  problem  treated  as  a  geometrical  transformation  on  the 
Poincar^  sphere. 

An  introduction  to  second  order  statistical  effects  in 
polarimetric  scattering  is  also  provided  via  the  coherency  matrix 
and  Stokes  vector.  The  matrix  governing  transformation  of  these 
second  order  parameters  is  related  to  the  elements  of  the 
coherent  scattering  matrix.  ^ 

This  memorandum  was  derived  from  a  set  of  lectures  given  by 
the  author  at  Birmingham  University  in  July  1983, 
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INTRODUCTION 


s 


When  an  object  of  arbitrary  size,  shape  and  material  scatters  an  incident 
electromagnetic  wave,  measurement  of  the  parameters  of  the  scattered  field  may 
be  used  to  yield  information  about  that  object.  One  very  important  parameter 
of  the  wave  which  has  not  been  fully  exploited  in  the  past  is  its  polarisation. 

On  the  whole  radars  have  measured  amplitude,  phase  and  doppler  shift  of  targets 
but  have  remained  fixed  polarisation  devices  and  thus  have  not  obtained  the 
maximum  available  information,  even  for  a  fixed  frequency  and  aspect.  The 
capability  of  obtaining  complete  polarimetric  information  is  therefore  an 
attractive  feature,  and  how  this  information  may  be  gathered  and  best  used  is 
the  subject  of  these  lectures. 

Although  polarimetric  techniques  have  found  widespread  use  in  optical 
instrumentation  and  materials  analysis  for  many  years,  the  relevance  of  polari¬ 
metric  phenomena  to  microwave  scattering  has  remained  primarily  a  topic  of 
theoretical  interest  only.  This  is  unfortunate  since  at  longer  wavelengths 
polarisation  is  a  more  significant  parameter  in  scattering  than  it  is  in  optics 
and  as  such  promises  to  yield  a  considerable  amount  of  target  information. 

Since  polarisation  may  be  described  as  the  spatially  directive  quality  of  a  wave 
then  knowledge  of  the  scattering  characteristics  of  an  object  might  be  expected 
to  yield  information  relating  to  target  symmetry,  geometry  and  material  properties. 
As  a  simple  example  consider  rain  clutter  suppression  radars  which  transmit  and 
receive  like  sense  circular  polarisation  on  the  basis  that  this  will  cancel  the 
backscatter  from  spherically  symmetrical  raindrops.  It  will  be  shown  in  the 
course  of  these  lectures  how  more  subtle  indicators  of  target  symmetry  may  be 
derived  from  knowledge  of  the  polarisation  transforming  properties  of  the  target. 

Recently  there  is  a  great  deal  of  interest  in  optimising  radar  system 
performance  and  in  target  identification  studies  both  of  which  require  a  more 
exact  target  descriptor  than  the  conventional  fixed  polarisation  radar  cross 
section^.  One  way  of  providing  this  is  to  dercribe  the  target  by  a  polarisation 
scattering  matrix  of  finite  dimensions  with  manipulation  of  this  matrix  in  the 
radar  signal  processing  yielding  target  information.  Essentially  this  means 
that  the  transmit/receive  antenna  assembly  no  longer  acts  as  a  polarisation 
filter  but  collects  full  polarimetric  information  and  allows  all  filtering  to  be 
performed  in  the  processing^.  With  ongoing  advances  in  the  field  of  digital 
and  analogue  signal  processing  technology  this  allows  for  practical  exploitation 
of  polarimetric  information. 
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The  importance  of  developing  a  full  polarimetric  theory  for  radar  scattering 
was  realised  in  the  late  nineteen  forties.  In  particular  E.M.  Kennaugh  demon¬ 
strated  the  importance  of  polarimetric  techniques  when,  in  a  series  of  reports 

(2) 

written  between  1949  and  1954  ,  he  introduced  the  important  concept  of  null  or 

(3) 

optimal  polarisations  for  radar  targets.  It  was  Sinclair  in  1948  who  was  the 

first  to  describe  the  polarisation  transformation  properties  by  a  2  x  2  coherent 

(4) 

scattering  matrix.  At  about  the  same  time  Mueller  developed  a  more  general 

matrix  calculus  for  handling  partially  polarised  waves  in  optics.  This  calculus, 

based  on  Stokes  vectors  was  applied  to  radar  scattering  much  later  by  J.R.  Huynen^ 

and  allows  a  polarimetric  theory  to  be  developed  for  handling  radar  scatter  from 

fluctuating  targets  like  chaff  or  clutter.  More  recent  developments  have  seen 

the  practical  implementation  of  programmable  polarisation  filter  design  and  its 

use  for  more  efficient  suppression  of  rain  clutter Significant  theoretical 

work  on  incorporating  polarimetric  techniques  with  broadband  scattering  theories 

(8) 

has  been  done  by  Boerner '  .  Despite  these  advances  the  full  potential  of  polari¬ 

metric  techniques  has  still  not  been  realised  and  one  can  only  speculate  as  to  the 
future.  Certainly  the  next  decade  will  see  the  emergence  of  more  and  more  radar 
systems  utilising  polarisation  for  optimised  performance. 

POLARISATION  DEFINITIONS 

Before  undertaking  a  full  polarimetric  description  of  target  scattering  it 
is  necessary  to  consider  the  various  methods  used  to  describe  the  polarisation 
of  a  wave.  In  this  chapter  we  shall  cover  three  methods:  firstly,  the  orthogonal 
component  decomposition  for  plane  monochromatic  waves.  The  geometrical  parameters 
of  the  polarisation  ellipse  will  then  be  related  to  these  plane  wave  parameters 
since  the  former  are  useful  when  considering  graphical  aids  to  polarisation 
problems.  Finally,  the  more  general  Stokes  vector  formalism  shall  be  developed 
and  related  to  the  other  representations  through  the  Poincar£  sphere. 

In  the  most  general  case  a  propagating  electromagnetic  wave  has  six  field 

components,  namely  three  mutually  perpendicular  electric  and  magnetic  field 

vectors  (E  E  E  H  H  H  ) .  For  most  radar  applications,  however,  it  is 
x,  y,  z,  x,  y,  z  rr  ’  ’ 

sufficient  to  consider  only  the  case  of  plane  time  harmonic  waves  travelling  in 
free  space  so  that  in  general  we  may  write  the  field  as 


'»,■  ,  v  lyy’jvv  ';*.~7v»r.  r.  ■*rrT.-TTT;,^.^T  ».* I.*  •'  -~*  ■**  •  *  •  ^  *  "'•l‘  *■*  "•*'•*■  »*  r-  •*-  •  •  *’*-'• 


E  "a  cos(ty  +  6)  -  la  E  ■  ri  H 

x  xx  x,y,  y,x, 


n  ■  impedance  of  free  space  ■  120  » 


a^  cos(\p  *  <5^)  -  lb  t|>  ■  wt  -  kr 


6,6  «  instantaneous  phase  of  E 

X  y  x,y 

components  at  t  *  0 


co 


angular  frequency 


.  2n 

k  "  T 


r  ■  range  from  source 

Here,  x,y,z  is  a  right  handed  set  of  cartesian  co-ordinates  with  the  direction 
of  propagation  of  the  wave  in  the  z  direction. 

Tf  the  plane  wave  assumption  is  relaxed  to  allow  partially  polarised  waves 
where  a  ,  a  ,  6  ,  6  are  functions  of  time  then  the  wave  must  be 

* y \  *  y  (9) 

specified  by  its  second  order  statistics  or  coherency  matrix 
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<E  E  *>  <E  E  *> 
xx  x  y 


<E  *E  >  <E  E  *> 
L  x  y  yy  ■ 


<>  *  time  average 


Note,  significantly,  that  J  is  Hermitian  and  so  may  be  written  in  the  general 
form 


[a  b-icT 

b+ic  d  J 

Thus  four  real  numbers  (a,b,c,d)  are  required  to  specify  the  second  order 
statistics  of  the  wave. 
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1.  ORTHOGONAL  COMPONENT  DECOMPOSITION 
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The  polarisation  of  a  plane  monochromatic  wave  is  defined  as  the  locus  of 
the  electric  field  vector  (Ex+Ey)  in  a  plane  perpendicular  to  the  direction  of 
propagation.  In  general  the  vector  direction  will  change  as  a  function  of  time, 
and  according  to  IEEE  standard  when  observed  along  a  line  of  sight  anti¬ 
parallel  to  the  direction  of  propagation  is  called  left  handed  polarisation  if 
it  rotates  clockwise  and  right  handed  if  counter  clockwise.  If  the  E  vector 
does  not  rotate  at  all,  the  wave  is  said  to  be  linearly  polarised. 


In  order  to  derive  an  expression  for  this  locus  for  arbitrary  a  ,a  ,6  ,<$ 

(9)  x  y  x  y 

combine  equations  la  and  lb  so  as  to  eliminate  il>  .  When  this  is  done  the 
equation  for  the  locus  is 


(S‘  •  0f  - 


E  E 

2  ■■  X  cos  <5 

Vy  ** 


•  2. 

sin  6 


6  -  6 
y  x 


This  is  the  equation  of  an  ellipse  for  which  linear  and  circular  polarisations  are 
special  cases.  In  general,  the  wave  may  be  represented  by  the  two  component 
column  vector 


E  j  x 
x 


E  :  y 
y 


xy 

a  e  :  y 

y  3 


J(t+6v) 


When  considering  only  the  state  of  polarisation  of  this  wave  then  use  will  be 
made  of  the  column  vector 


Formally,  this  vector,  each  element  of  which  is  complex,  forms  a  spinor  and  much 
of  what  follows  finds  elegant  interpretation  in  terms  of  spinor  algebra^11^. 

Note  that  in  the  above  notation  x  and  y  may  be  any  two  orthogonal  polarisations 
of  which  horizontal  and  vertical  and  left  and  right  circular  are  common  examples. 
In  general  any  pair  of  elliptical  polarisations  will  suffice  as  long  as  they 


4 


satisfy  the  orthogonality  condition 


x  •  y  "  0 

Each  pair  forms  a  so-called  polarisation  base  written  (x,y).  Any  polarisation  P 
may  then  be  represented  in  terms  of  a  complex  combination  of  these  components 

x,y  -  unit  orthogonal  vectors 
A,B  -  complex  coefficients  determining  P 

Of  interest  is  to  consider  how  to  express  P  in  terms  of  an  arbitrary  base  (x',y*) 
when  it  is  known  in  terms  of  (x,y).  The  most  general  transformation  (x',y')  from 
(x,y)  has  the  following  normalisation  and  orthogonality  conditions 

x'  •  x’*  ■  1 
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This  leads  to  the  conditions 


$1  “  *3  "  *2  “  ^4 


e 


Further,  we  can  set  ■  0  without  loss  of  generality  and  so  can  write  the 

transformation  in  matrix  notation  as 


^  l" 

r  .  J6  -i 

“  A- 

x1 

SR 

cosa  sma  e 

X 

*  9 

-J<5 

A 

LyJ 

-sina  e  cosa 

_y_ 

This  may  also  be  written 


[T] 


cosa 


1  tana  e 

-J<5 


J6 


-tana  e 
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At  this  point  is  is  useful  to  introduce  the  complex  polarisation  ratio  defined 
as 

B  *14* 

P  "  A  ■  M  e  p 

This  ratio  defines  P  in  terms  of  the  ratio  of  the  complex  coefficients  A  and  B. 
From  this  definition  and  the  general  form  of  [  T]  it  is  evident  that 


p  ■  tana  e 


J6 


where 


tana 


0  <  a  <  90 


*B  '  *A 


0  <  6  <  360 


The  geometrical  significance  of  the  two  angles  a, 6  will  be  discussed  later  in 
this  chapter  when  considering  the  Poincar£  sphere.  Note  that  in  terms  of  p 


cosa  ■  (1  +  pp*) 


-J 
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This  then  allows  us  to  write  [T]  ,  the  transpose  of  [T] ,  which  governs  the  trans¬ 

formation  of  coefficients  of  (x  y)  as 


1  -p* 


P  1 


2.  GEOMETRICAL  PARAMETERS 


A  convenient  method  of  representation  of  the  polarisation  ellipse  is  in 
terms  of  its  inclination  angle  0  (0  <  0  <  180°)  and  ellipticity  angle  t 
(-45°  <  t  <  45°)  as  defined  in  Figure  1.  By  convention,  positive  values  of  t 
correspond  to  left  hand  polarisations  and  negative  values  right  hand.  The 
amplitude  of  the  wave  is  defined  as  in  Figure  1. 

In  order  to  express  an  arbitrary  state  of  polarisation  P  in  matrix  form 
consider  firstly  that  the  x,y  axes  in  Figure  1  are  aligned  with  the  major  and 
minor  axes  of  the  ellipse  respectively.  From  the  definition  of  x  we  can  write 


a  cos  ut 
x 


a  cos(ut+5  ) 

y  xy 


a  cost  coswt 


-a  sinx  sinwt 


a  cost 

Re  e“ 

Ja  sinx 


In  complex  notation 


[cost  "1 
J  sintj 


P  -  a 


For  the  general  case  where  the  major  axis  of  the  ellipse  makes  an  angle  0  with 
the  axes,  this  column  vector  is  multiplied  by  the  transformation  matrix  for 
rotations  in  2-dimensions  to  yield 


:os6  -s 


P  -  a 


iin0  cos 


in©"jr  cost  "J 
OS0JLJ  sinxj 


P(a,0,x) 


In  terms  of  these  geometrical  parameters,  the  orthogonal  polarisation  to  P(a,0,T) 
is  defined  as 
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Pj_  "  P(a,  e  +  j,  -t) 


In  other  words  an  ellipse  with  its  major  axis  rotated  through  90  and  with  the 
same  magnitude  of  ellipticity  but  opposite  sense  (Figure  2).  The  transformation 
matrix  [  Tj  changing  (x,y)  into  (x’.y')  comprises  two  components:  one  representing 
a  deformation  of  ellipticity  given  by  angle  and  the  other  a  rotation  of  the 
major  axis  by  0c.  Thus 


[T]  - 


|"cos0  -sin0 
c  c 


sin0  cos© 

C  C‘ 


cost  J  sini 
c  c 


J  sini  cost 

c  c 


By  expanding  this  matrix  product  and  comparing  terms  with  [  T]  written  in  terms 
of  a ,6,  the  relationships  between  the  geometrical  parameters  and  the  polarisation 
ratio  terms  may  be  derived  as^^ 


cos2a 


cos20  cos2t 


tan6  «  tan2t  csc20 


with  the  inverse  relationships 


tan20  *  tan2ot  cos6 


a 

R 
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sin2t  *  sin2a  sin6 


STOKES  PARAMETERS  AND  THE  COHERENCY  MATRIX 


The  above  definitions  are  adequate  when  considering  coherent  waves 
ie  waves  where  the  parameters  of  the  polarisation  ellipse  are  independent  of  time. 
More  generally  it  is  of  interest  to  consider  scattering  from  fluctuating  targets 
and  in  these  instances  the  polarisation  ellipse  will  fluctuate  as  a  function  of 

time.  In  particular  <E  >  and  <E  >  may  equal  zero  for  uniformly  random  variations 

.  .  x  y  .  (9) 

m  phase  and  so  in  order  to  describe  the  wave  one  must  use  its  coherency  matrix 

This  measures  the  complex  correlation  that  exists  between  two  orthogonal  polari¬ 
sation  vectors  over  the  period  of  observation.  At  the  two  extremes  there  will  be 
either  complete  correlation  in  which  case  the  wave  is  coherent  and  elliptically 
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polarised  (ep)  or  zero  correlation  which  implies  a  randomly  polarised  wave  (rp). 
In  general  there  will  exist  some  degree  of  correlation  and  hence  the  wave  is 
termed  partially  polarised  (pp) . 


Previously  the  coherency  matrix  was  defined  as 


<E  E  *> 

X  x 

<E  E  *> 

x  y 

<a 

2> 

X 

-J<5 

<a  a  e 

x  y 

<E  *E  > 

L  x  y 

<E  E  *> 

y  y  J 

<a  a 

L  x  y 

J6 

e  xy> 

<a  2> 

y  -J 

a 

b-ic 

"j 

XX 

J 

xy 

b+ic 

d 

J 

-  yx 

j 

yyJ 

Note  that: 


Trace  (J)  -  sum  of  diagonal  elements  -  Jxx  +  -  total  power  in  the  wave. 


Also,  since  J  J  are  non-negative: 
xx  yy 


det  (J)  *  J  J  -  J  J  >0 
xx  yy  xy  yx 


For  rp  waves  there  are  no  preferred  polarisation  parameters  and 


so 


"  J  and  J  -  J  «0 
xx  yy  xy  yx  v 


[  J]  -  I 


t: :] 


At  the  other  extreme,  for  ep  waves  a  a  6  are  time  independent  so 

x  y  xy 


det  (J)  -  a  a  2  -  a  2a  2  e  xy  e  ** 
x  y  x  y 


(  J] 


■r  i 

Lq*  rJ 


PR  -  QQ*  -  0 
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Of  particular  interest  is  the  decomposition  of  a  partially  polarised  wave  into 
the  sum  of  an  rp  and  and  ep  wave 


J  -  J  +  J 
pp  rp  ep 


where 


J 

”Jxx 

J  1 

xy 

j  ■ 

'i 

0 

j  * 

‘p  q' 

pp 

J 

J 

rp 

0 

i 

ep 

Q*  R 

L  yx 

yyJ 

U  -J 

L  J 

P,Q,R,I  >  0 


PR  -  QQ 


xx 


I  +  P 


J  -  Q 
xy 


J  -  Q 
yx 


J  -  I  +  R 

yy 


From  the  determinant  equation 


(J  -I)(J  -I)  -  J  J  -  0 

xx  yy  xy  yx 


Thus  I  is  an  eigenvalue  of  the  coherency  matrix  given  by 


I  -  i(J  +  J  )  t  HU  +  J  )2  -  4  det  (J))* 
*  xx  yy  *  xx  yy 


Both  roots  of  this  equation  are  real  and  non-negative  but  since  P,R  >  0  the 
negative  sign  must  be  taken  showing  that  the  decomposition  is  unique.  As  before 
the  total  intensity  of  the  wave  is  given  by  Tr(J)  ■  J*x  +  J  and  the  total  power 
density  in  the  ep  part  of  the  decomposition 

T  (J  )  -  P  +  R  -  ((J  +  J  )2  -  4  det(J))* 

r  ep  ,:x  yy 

From  this  the  degree  of  polarisation  of  the  wave  is  defined  as 

D  «  power  in  ep  component 
p  total  power 


1  - 


4  det (J) 

(J  +J  )‘ 
xx  yy 


STOKES  PARAMETERS 


Of  particular  interest  is  the  representation  of  the  coherency  properties  by 

four  real  quantities  instead  of  by  the  complex  correlation  matrix.  One  such 

choice  of  quantities  was  developed  by  Sir  George  Stokes  in  1852  for  use  in 

vibrational  theories  of  light  propagation  in  the  ether.  They  are  formally 
(13) 

defined  as  the  four  components  of  the  associated  longitudinal  vector  of 
the  polarisation  spinor  P  as  defined  previously  and  as  such  allow  the  use  of  the 

(14) 

Minkowski  model  of  Lorentz  space  for  handling  polarisation  problems.  They 
are  related  to  the  elements  of  the  coherency  matrix  as  follows: 


80 

-  J  (J 
*  xx 

+  J  )  - 

yy 

J(a  +  d) 

-  T  (J) 
r 

*  <a 

81 

-  |  (J 

2  XX 

-  J  )  - 

yy 

i (a  -  d> 

-  <ax2> 

-  <a?2 

82 

-  J(J 

2  xy 

+  J  )  n 
yx 

b 

“  <Sla2 

cosfi^ 

g3 

*  ii  (J 

5  v  xy 

-  J  )  ■ 

yx7 

c 

*  <ala2 

sindj^ 

(gQ  g2  83)  ^  a  Stokes  vector 


Physical  interpretation  may  be  placed  on  the  elements  of  £  by  considering  the 
case  where  the  coherency  matrix  is  represented  in  terms  of  (h,v)  base.  Under 
these  circumstances  g^  represents  the  total  power  density  of  the  wave;  is  a 
measure  of  how  much  like  vertical  or  horizontal  the  wave  polarisation  is;  gj  a 
similar  measure  of  ±  45°  nature  and  a  measure  of  the  ellipticity  of  the 
polarisation. 


A  derivation  of  the  Stokes  parameters  from  the  coherency  matrix  may  be  had 

(15) 

by  considering  an  expansion  of  the  latter  in  terms  of  the  Pauli  spin  matrices 
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These  four  matrices  form  a  representation  of  a  mathematical  group  with  the 
following  multiplication  table: 
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It  is  well  known that  any  2x2  matrix  [  S]  may  be  expanded  in  terms  of  this 
set  so  that 
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where  are  the  coefficients  of  expansion  given  by 
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For  the  particular  case  of  the  coherency  matrix  the  expansion  yields 
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Notice  the  change  of  order  of  Stokes  parameters  when  compared  with  the  previous 
definition,  an  unfortunate  complication  arising  in  many  areas  of  polarisation 
algebra  and  one  which  causes  untold  confusion  when  carried  through  to  scattering 
matrix  theory.  The  above  is  known  as  the  natural  ordering  of  Stokes  parameters 
and  is  frequently  used  in  optics  whereas  the  previous  definition  is  termed 
traditional.  Unfortunately  in  the  radar  literature^  there  is  yet  another 
permutation  used  and  for  want  of  a  better  name  I  shall  call  this  the  radar 
ordering.  This  arises  from  defining  the  Pauli  matrices  as 


3 


k.i 

n 


V 

-A 


Notice  that  in  this  case  S  must  be  multiplied  by 


0  -i 


in  order  to  achieve  real  Stokes  parameters. 


In  the  radar  ordering  the  Stokes  vector  takes  the  form 


£  m  ^8g*  83 »  82*  8j) 


To  transform  between  the  different  systems  the  following  matrices  must  be  used 
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Occasionally  the  so-called  modified  Stokes  vectors  are  used  which  are  related 
to  £  -  (80.81»82»g3)  by 
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The  orthogonality  condition  for  Stokes  vectors  and  g  is 
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Only  the  traditional  ordering  shall  be  considered  in  the  remainder  of  this  chapter 
since  it  leads  most  conveniently  to  the  Poincard  sphere  representation  of 
polarisation. 

Comparison  of  the  properties  of  the  Stokes  vector  with  the  coherency  matrix 
gives  rise  to  the  following  useful  relationships. 

The  determinant  condition  for  J  becomes 
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This  is  known  as  the  condition  for  physical  realisability  of  the  Stokes  vector 
and  for  an  ep  wave  becomes 
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For  an  rp  wave 
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The  decomposition  theorem  becomes 
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Finally,  the  degree  of  polarisation  is  given  by 
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The  condition  for  physical  realisability  states  that  for  an  ep  wave 
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The  three  Stokes  parameters  (g^  %2  S3)  may  be  thought  of  as  the  co-ordinates  of 
a  point  in  three  dimensional  space  with  distance  from  the  origin  given  by  gQ. 
Further,  this  radius  is  equal  to  the  amplitude  of  the  wave  and  as  such  the  loci 
of  polarisations  of  equal  amplitude  is  a  sphere.  In  optics  this  is  called  the 
Poinar£  sphere,  there  being  a  one  to  one  correspondence  between  points  on  the 
surface  and  the  set  of  all  possible  polarisations.  Note  that  this  spherical 
geometry  is  only  true  for  coherent  scattering  and  for  partially  polarised  waves 
the  transformations  occur  on  different  spheres  .  it  is  a  feature  of  the 
Stokes  vector  formalism  that  it  allows  the  treatment  of  such  problems  and  as 
such  forms  the  basis  for  a  general  theory  of  depolarisation  for  which  coherent 
scattering  is  a  special  case.  Shown  in  Figure  3  are  some  examples  of  normalised 
Stokes  vectors  for  some  commonly  used  polarisations. 


The  transformation  equations  governing  the  movement  of  a  point  over  the 
surface  of  a  sphere  are  well  known  from  trigonometry^17^  and  are  best  represented 
by  three  Euler  angles,  each  representing  a  rotation  about  the  x,  y,  z  axis 
respectively.  The  geometry  of  this  problem  is  shown  in  Figure  3  which  shows  a 
right  handed  cartesian  co-ordinate  system  together  with  various  co-ordinates  of 

A  A 

the  point  P  in  polarisation  space.  For  transformation  of  P  by  an  angle  v  about 
the  x-axis  the  governing  transformation  matrix  for  £  is 
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Similarly  for  transformations  by  t  and  \p  about  the  y  and  z  axes  respectively 
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A  general  transformation  is  then  written 


R  -  R2(*)  Ry (t)  Rx(v) 


Rather  than  consider  the  cartesian  co-ordinates  of  P  it  is  sometimes  convenient 
to  consider  the  angular  co-ordinates  (latitude  And  longitude).  In  order  to  see 
how  these  are  related  to  the  geometrical  parameters  of  the  polarisation  ellipse 
consider  again  the  original  form  of  the  coherency  matrix. 


<E  E  *>  <E  E  *> 
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In  terms  of  the  geometrical  parameters 


cose  cost  -  isin6  sint 


sin6  cost  +  icos6  sinT 


<E  E  *> 
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<E  E  *> 
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g^  ■  a  cos20  cos2x 

2 

g2  ■  a  sin26  cos2x 

2 

g^  ■  a  sin2x 

From  standard  spherical  trigonometry  these  are  seen  to  be  the  equations  relating 
spherical  polar  to  cartesian  co-ordinates  and  hence  the  remarkable  result  that 
the  latitude  and  longitude  of  a  point  P  on  the  Poincar€  sphere  are  2x  and  26 
respectively  where  6  and  x  are  the  inclination  angle  and  ellipticity  angle  of 
the  corresponding  polarisation  ellipse  (Figure  4). 

The  relations  between  0,x  and  a, 6  (the  polarisation  ratio  parameters)  are 

(18) 

again  standard  results  from  spherical  trigonometry v  .  They  form  the  elements 
of  a  spherical  triangle  as  shown  in  Figure  5.  The  a, 6  co-ordinates  figure 
prominently  in  the  theory  of  null  polarisations. 

The  Poincard  sphere  has  the  following  interesting  properties: 

a.  The  poles  of  the  sphere  represent  left  and  right  circular  polarisations 
(x  •  ±  45°) . 

b.  The  upper  and  lower  hemispheres  map  similar  sets  of  elliptical  polari¬ 
sations  with  opposite  sense.  The  upper  hemisphere  is  chosen  arbitrarily  for 
left  sense  polarisations  (positive  x) . 

c.  The  loci  of  polarisations  of  equal  ellipticity  are  in  planes  parallel 
to  the  equator  which  itself  represents  the  set  of  all  linear  polarisations. 

d.  Orthogonal  polarisations  lie  diametrically  opposite  on  the  sphere. 

-It  is  worth  noting  the  similarity  between  the  Poincar£  sphere  and  its  plane 
(19) 

projections  and  a  similar  geometry  used  for  impedance  calculation  in  circuit 
theory.  Indeed  the  Smith  chart  is  a  projection  of  the  impedance  sphere  and  may 
be  similarly  used  to  solve  polarisation  problems.  This  commonality  between  two 
seemingly  different  areas  of  work  can  be  quite  useful  for  understanding 
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(19) 

transformations  on  the  sphere  .  The  similarity  arises  out  of  the  fact  that 
p  is  a  complex  ratio  of  like  field  components  in  a  plane  perpendicular  to 
direction  of  propagation  while  the  impedance  Z  is  a  similar  ratio  of  electric 
and  magnetic  field  components. 

Although  the  Smith  chart  could  be  used  for  polarisation  problems  the  most 
common  plane  projection  used  is  the  Polarisation  Chart  as  shown  in  Figure  6. 

Note  that  two  such  charts  are  required  to  allow  for  mappings  on  both  hemispheres, 
the  chart  being  a  projection  of  the  sphere  onto  the  equatorial  plane  so  that  linear 
polarisations  lie  around  the  circumference  with  circular  in  the  centre.  There 
has  also  been  a  scale  change  in  order  to  give  a  linear  reduction  in  eccentricity 
from  unity  in  the  centre  to  zero  on  the  circumference. 

RECEPTION  OF  POLARISED  WAVES  BY  AN  ANTENNA 

In  this  section  an  expression  will  be  derived  for  the  power  received  by  an 

antenna  whose  polarisation  is  fixed  as  P,  being  used  to  receive  a  plane  wave  with 
* 

polarisation  Q.  An  ideal  antenna  will  be  assumed  for  the  sake  of  clarity  and  the 
more  general  problem  of  non-ideal  antenna  characteristics  discussed  in  a  later 
chapter. 

The  polarisation  of  an  antenna  is  defined  as  the  polarisation  of  the  wave 
it  radiates  on  transmit  (remembering  that  this  is  defined  in  a  right  hand  set  of 
x,y,z  co-ordinates  where  the  positive  z  direction  is  pointing  in  the  direction 
of  propagation  of  the  wave).  Using  the  geometrical  parameter  representation  we 
may  then  write  the  antenna  polarisation  as 


“cose 

-sine 

COST 

A 

p 

P 

P 

P  -  a 

P  P 

sine 

D 

cose 

o 

isinr 

D 

’a  '  is  the  gain  function  of  the  antenna.  The  geometry  of  the  problem  to  be 

.  A 

solved  is  then  as  shown  in  Figure  7.  The  important  point  to  note  is  that  Q  is 
defined  in  its  own  set  of  right  handed  co-ordinates  which  are  different  to  those 
of  P.  One  system  is  obtained  from  the  other  by  a  180  rotation  about  the  y  axis. 
In  matrix  notation  the  x,y  co-ordinates  may  be  related  as 
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Although  it  is  desirable  to  define  the  polarisation  of  an  incoming  wave  in  terms 
of  its  own  co-ordinate  system  it  is  also  convenient  to  derive  an  expression  for 
the  received  power  using  polarisations  defined  in  the  same  co-ordinates  since 
then  the  same  point  on  the  Poincard  sphere  may  be  used  to  represent  an  antenna 
whether  it  transmits  or  receives. 


In  its  own  co-ordinate  system  the  incoming  wave  is  defined  by 
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In  order  to  express  this  in  terms  of  the  receiving  antenna's  co-ordinates 
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It  is  easily  verified  that 


*  Qn(aa»"0n»  "O 
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Thus,  if  a  wave  of  polarisation  is  incident  upon  an  antenna,  then  for  that 
antenna  to  be  matched  to  the  wave  in  the  x,y  plane  it  must  not  have  the  same 
polarisation  as  but  the  one  given  as  ,  ie  opposite  sense  and  negative 
inclination  angle.  This  is  called  the  symmetric  polarisation  to  and  is  shown 
geometrically  in  Figure  8.  This  may  seem  a  remarkable  theorem  and  has  certainly 
been  at  the  root  of  much  confusion  in  the  past  but  is  easily  verified  by 
considering  transmission  between  two  facing  antennae,  both  of  which  are  polarised 
at  45°  linear.  In  this  case  they  are  orthogonally  polarised  to  each  other  and  so 
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no  power  will  be  received  at  all.  In  fact  it  is  only  for  vertical  and  horizontal 
that  the  symmetric  polarisation  is  equal  to  its  parent. 
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Also  shown  in  Figure  8  are  two  other  useful  polarisations  related  to  the 

A 

parent  Q.  The  conji 
of  sense  and  is  given  by 


same  parent  Q.  The  conjugate  polarisation  is  identical  to  except  for  a  change 


6  (a  ,8  ) 

c  c9  c 


The  transverse  polarisation  has  the  same  sense  but  negative  the  inclination  angle 

VWV 

These  polarisations  are  useful  when  considering  the  backscatter  from  targets  with 
certain  symmetries.  Figure  9  shows  how  these  polarisations  relate  to  their 
parent  on  the  Poincar£  sphere.  It  is  convenient  that  by  adopting  the  co-ordinate 
change  described  we  can  dispense  with  such  transformations on  the  sphere  and  write 
the  equation  for  received  voltage  as 


V  «  P  •  Q  -  p  .  Q 
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where  P  and  Q  are  defined  in  their  own  co-ordinate  systems  and  the  dot  product 
for  complex  vectors  is  defined  as 
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Note  that  from  this  definition 
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Before  working  out  an  explicit  form  for  V  it  is  worth  noting  a  useful  form  of 
polarisation  matrix  algebra,  based  on  the  Pauli  matrices  discussed  earlier, 
which  allow  a  short  hand  approach  to  be  used  in  solving  what  would  otherwise 
lengthy  trigonometrical  calculations.  We  know 


The  6  matrix  may  be  expanded  in  terms  of  the  Pauli  matrices  to  yield 
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Similarly  for  the  t  dependence 
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And  for  the  sake  of  completeness 
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This  extension  of  the  exponential  notation  from  ordinary  complex  number  theory 
to  Pauli  matrices  is  a  powerful  technique  and  allows  the  general  polarisation 
ellipse  to  be  written  as 
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The  following  rules  are  useful  for  exploiting  this  algebra: 


Proof 


(a+0)c. 


cos(ci+$)Oq  +  isin(a+g)o2 


ao0  $o0 


e  e  *■  -  (coscxOq  +  isinac^)  (cosBct^  +  ising02) 


cosa  cosB  +  icosa  sini^  +  isina  COS6O2  ~  sina  sinB 


cos(a+6)  +  isin(a+B)o, 


a a,  -00 

2.  Oj  e  *  -  e  1  0j 


Proof 


o^(cosaOQ  +  isinaa2)  ■  (cosao^  -  isinao2)o^  (see  multiplication  table) 
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Using  this  notation  the  expression  for  received  voltage 
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And  considering  only  the  central  portion 
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In  matrix  form  this  becomes 
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and  when  multiplied  out  the  expression  for  V  becomes 
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The  received  power  is  given  by  W* 
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and  by  using  the  relationships 
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This  is  a  fundamental  equation  in  polarisation  analyses  and  as  such  warrants 
further  comment.  Note  that  for  maximum  power  received  the  conditions 


0  *  -6  and  t  ■  t  must  be  satisfied. 

p  q  p  q 


In  this  case 
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Notice  that  the  condition  on  the  ellipticity  is  not  as  might  be  expected  since 
it  states  that  maximum  power  is  received  when  the  incident  wave  has  the  same 
sense  as  the  antenna,  when  both  are  quoted  in  their  own  co-ordinate  system.  This 
is  due  to  the  fact  that  the  sense  of  polarisation  changes  not  only  with  the 
rotation  of  co-ordinates  provided  by  but  also  depending  on  the  direction  of 
observation,  ie  either  parallel  or  antiparallel  to  the  direction  of  propagation. 


There  are  four  independent  power  measurements  required  by  an  antenna  in 
order  to  determine  Q  as  evidenced  by  the  equation  for  P.  Usually  the  four  used 
are  linear  vertical  and  horizontal,  '-A5°  linear  and  a  circular  polarisation^21^ 
and  this  allows  the  determination  of  (6^  t^) .  Alternatively  a  dual  channel 
receiver  could  be  used  measuring  simultaneously  the  vertical  and  horizontal 
components  of  Q  and  as  long  as  some  measure  of.  the  time  phase  angle  between  V 
and  H  is  provided  this  allows  for  the  instantaneous  measurement  of  polarisation. 
This  latter  technique  will  be  used  later  in  order  to  determine  the  scattering 
from  targets. 

In  terms  of  the  Stokes  parameters  the  equation  for  P  takes  on  the  form  of 
a  dot  product,  namely 


where 


P  *  80h0  +  8lhl  ’  g2h2  +  g3h3 


£  *  (8q  82  83)  “  Stokes  vector  for  Q,  the  incident  wave 

h  *  (h^  h^  h 2  h^)  “  Stokes  vector  for  P,  the  antenna  polarisation 
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Notice  that  this  may  be  written  in  matrix  form  as 


P  - 


1  0 
0  1 
0  0 
0  0 


0 

0 

-1 

0 


0 

0 

0 

1 


IM]  & 


g' 


In  general  the  received  power  P  will  not  be  a  maximum  and  this  leads  to  the 


definition  of  the  polarisation  efficiency 


r  - 


2  2 
a  a 
P  q 


f  ■  1  for  polarisation  matched  antenna 


T  ■  0  for  orthogonally  polarised  antenna 


A  A 

On  the  Poincarg  sphere  P  may  be  mapped  as  the  antenna  polarisation  and  then  P^ 
lies  diametrically  opposite.  If  Q  is  the  polarisation  of  the  incident  wave  the 
polarisation  efficiency  may  bo  related  quite  simply  to  the  angular  separation 
between  P  and  Q  as  follows.  Figure  10  shows  the  two  polarisations  mapped  on  the 
sphere  with  2 the  angular  separation.  Figure  11  shows  the  so-called  power 
density  semicircle  from  which  it  is  evident  that  the  normalised  power  in  Q  is 


2  2 
cos  4>  *  sin  i it 


A  A  A  A 

If  Q  is  then  decomposed  into  (P  P,  ),  the  fractional  power  contained  in  the  P  and 

2  2  A 

components  is  cos  v  and  sin  ij;  respectively.  From  the  definition  of  T  this 
means 


2, 

cos  ^ 


This  relationship  is  of  great  interest  in  polarisation  filter  design  since  it 
clearly  demonstrates  that  the  loci  of  constant  response  are  circles  centred  on 
the  antenna  polarisation.  One  of  the  reasons  the  polarisation  chart  is  used  as 
a  plane  projection  is  that  circles  on  the  sphere  translate  to  circles  on  the 


chart  and  so  allow  the  design  of  filters  with  the  desired  response  in  polarisation 
(7) 

space  . 
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TARGET  SCATTERING  MATRIX 


In  general  terms  a  radar  target  may  be  described  as  a  polarisation  trans¬ 
former  in  that  it  operates  on  the  incident  wave  such  that  the  polarisation  of 
the  reflected  wave  bears  a  complicated  but  deterministic  relationship  to  the 
incident  polarisation.  Further,  the  nature  of  this  transformation  will  depend 
on  the  geometry,  surface  structure  and  material  composition  of  the  target  and  in 
the  remaining  part  of  these  lectures  it  will  be  shown  how  by  making  appropriate 
field  measurements  and  processing  the  data  in  the  correct  way  this  information 
may  be  made  available  for  enhanced  target  detection  and  target  identification. 

For  a  complex  radar  target  like  an  aircraft  the  scattering  will  not  only  be 
a  function  of  incident  polarisation  but  of  target  aspect  and  radar  frequency. 

It  will  be  assumed  initially  that  these  are  fixed  and  that  the  scattered  field 
parameters  are  linearly  related  to  the  incident  field.  The  polarisation  trans¬ 
forming  properties  may  then  be  expressed  as  a  transformation  matrix,  the  polari¬ 
sation  scattering  matrix.  The  dimensions  and  form  of  this  matrix  depend  on  the 
representation  cf  polarisation  used:  for  coherent  scattering  the  incident  and 
reflected  waves  may  be  expressed  in  terms  of  polarisation  base  (x,y)  and  the 
matrix  is  then  2x2  complex  (Figure  12).  Note  that  in  general  eight  measure¬ 
ments  will  be  required;  four  amplitudes  and  four  phase  angles,  for  the  determina¬ 
tion  of  this  matrix. 

If  the  Stokes  parameters  had  been  used  then  the  scattering  matrix  would  be 
4x4  and  real  (Figure  12).  This  then  requires  sixteen  amplitude  measurements 
to  be  made  but  will  contain  extra  information  above  that  provided  by  the  coherent 
matrix  since  it  is  based  on  the  more  general  Stokes  vector  formalism.  This 
matrix  is  known  as  the  Stokes  reflection  or  Mueller  matrix  and  is  of  fundamental 
importance  in  describing  the  scattering  of  partially  polarised  waves.  In  1970, 
Huynen^  ^  developed  a  decomposition  theorem  analagous  to  the  decomposition  of  the 
coherency  matrix  enabling  a  set  of  average  coherent  scattering  matrix  parameters 
to  be  obtained  from  the  Stokes  matrix.  In  this  chapter  we  shall  consider  only 
the  coherent  matri::  and  its  properties. 


S'  W  WW.W 
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transmit 


‘  ISJ  *t 


It  is  assumed  that  [  S]  in  this  form  has  been  premultiplied  by  a ^  so  that  the  h,v 
axes  are  the  same  for  transmit  and  receive  as  required  by  the  Poincarg  sphere 

* 

representation.  In  the  notation  employed  S  means  transmitting  polarisation  x 
and  receiving  polarisation  y  in  both  amplitude  and  phase. 


The  voltage  received  at  a  pair  of  antenna  terminals  for  fixed  transmit 
polarisation  P^t  receiver  polarisation  PR  and  target  scattering  matrix  [  S]  may 
be  written 


V  -  #R  *  l  S]  pt 

For  monostatic  radar  systems  (ie  ones  using  the  same  antenna  for  transmit 
and  receive)  various  simplifications  may  be  made  to  [  S]  . 

The  absolute  phase  of  the  target  is  a  function  of  its  range  and  velocity 
and  as  such  is  not  a  target  related  parameter  (unless  doppler  type  analysis  is 
required).  Therefore  there  is  no  loss  of  generality  if  one  of  the  phase  angles 
in  [  S]  is  set  to  zero  and  the  others  measured  relative  to  the  corresponding 
phase  centre  on  the  target.  This  then  reduces  the  number  of  required  measurables 
to  seven. 

For  monostatic  systems  the  reciprocity  theorem  for  antennas  demands 

V  -  PA  •  [SJ  PB  -  Pg  •  [S]  PA 

but 

PA  •  IS]  PB  -  is]1  PA  •  PB  -  PB  •  rsjT  PA 
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which  implies 

(S]T  -  [S] 

This  means  that  for  monostatic  systems  the  scattering  matrix  is  symmetric 


le 


HV 


-  S 


VH 


If  the  phase  of  these  diagonal  terms  is  taken  as  reference  then  the  so-called 
relative  scattering  matrix  results 


V 

•|sHHl  lsHV! 
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receive 
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_  v_ 

transmit 

This  matrix  is  determined  by  measurement  of  five  parameters,  namely  three 
amplitudes  and  two  relative  phase  angles. 

Occasionally  another  2x2  target  matrix  is  used,  called  Grave's  Fower 
matrix  it  is  related  to  [  S]  as  follows: 

£s  -  ( s]  It 

The  scattered  power  is  given  by 

£  *  *  Ig  *  £  *T  I  S]  *T[  S]  £t  -  E*T  P  £t 
"a  c" 

P  -  S*TS  * 

_c*  b. 

T* 

where  S  is  the  transpose  conjugate  of  [  S]  and  P  is  Hermitian.  For  example,  in 
(h,v) 
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C  ■  SHH  SHV  +  SVH  SW 


Thus  's'  is  the  total  available  backscattered  power  in  the  horizontal  component 
due  to  both  the  SHH  and  S^y  elements  of  {  S] .  Measurement  of  this  matrix  may  be 
had  by  transmitting  four  different  polarisations  and  receiving  the  total  power 
backscattered  in  each  case.  This  matrix  contains  a  subset  of  the  target  informa¬ 
tion  contained  in  [  S]  and  so  will  not  be  considered  further  in  these  lectures. 

Figure  13  shows  some  examples  of  the  relative  phase  scattering  matrix  for 
simple  radar  targets.  Note  that  all  these  matrices  are  quoted  in  (h,v). 

The  identity  matrix  is  indicative  of  reflection  from  a  flat  plate  at  normal 
incidence  or  from  a  trihedral  retroref ] ector  (both  are  same  in  a  polarimetric 
sense)  while  the  double  bounce  dihedral  reflector  has  a  180°  phase  difference 
between  HH  and  W.  Also  shown  are  some  other  important  scattering  types  such 
as  the  linear  target  and  helix  both  of  which  are  important  when  considering  the 
classification  of  targets  as  will  be  covered  later  in  this  course. 

The  most  powerful  aspect  of  scattering  matrix  measurement  is  that  it  allows 
the  prediction  of  target  scattering  for  any  transmit  polarisation  and  as  such 
provides  full  polarimetric  information  about  the  target. 

The  transformation  equations  governing  the  prediction  of  target  backscatter 
in  base  (x,y)  when  it  is  known  in  (x'.y*)  may  be  developed  using  the  change  of 
base  matrix  [  T]  derived  earlier.  Care  roust  be  taken  however  to  remember  that 
[  S]  is  an  operator  relating  an  incident  to  a  backscattered  wave  and  as  such, 
co-ordinate  changes  must  be  taken  into  account.  The  following  relationships 
hold: 


1.  tT]  is  a  unitary  matrix  operator  with  unit  determinant 


det(T)  -  1 
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2.  For  the  incident  system 

£t(x,y)  -  [T]  Zt(x', yf) 

3.  For  the  scattered  system 

£s(x,y)  -  IT]*  ^(x'.y1) 

The  conjugate  operator  is  employed  because  there  is  a  change  of  sense  of 
polarisation  when  the  scattered  co-ordinate  system  is  observed  from  the  trans¬ 
mitter  co-ordinates. 

From  the  definition  of  the  scattering  matrix 
Es(x,y)  -  IS(x,y)]  £t(x,y) 

-  [S(x,y)][T]  Et(x*,y') 

From  3 

£s(x\y’)  -  [T]T  £s(x',y') 

Is(x'  ,y')  -  [T]T(S(x,y)][T]  £t(x,y) 

Thus  in  general  the  transformation  of  [  S]  is  given  by  the  congruent  transformation 
[S(x’,y’)]  -  [T]T[S<x,y)][T] 


By  putting 


IS] 
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and  expanding  the  matrix  product,  the  transformation  equations  become 


SX*X'  "  (1+pp*>  (p  SYY  4  p(SXY  +  SYXJ  +  SXXJ 


SX'Y ' 


(1+PP*)  1  [p  -  P*5}^  +  S  yry  “  PP^S^^l 


SY*X’ 


C1*PP*)-1  Ip  Syy  -  p*s^x  *  svx  -  pp's^l 


SY’Y' 


U+pp*)'1  I  »%,  ♦  -  p'lS^  <  SYX)] 


For  monostatic  scattering  these  simplify  to 


SX'X’  "  ^1+pp  ^  ^ p  SYY  +  2SXYP  +  SXX^ 


Sx,y,  -  (1+ppV1  IpSyy  -  p*Sxx  +  S^U-PP*)] 
SY'X'  "  (1+PP*)  1  1  PSyy  ”  P*SXX  4  Sjjyd-PP*)] 


Sy'Y* 


(1+PP*)  1  [p  2SXX  +  SyY  -  P*  2S^y] 


Note  that  the  symmetry  of  the  monostatic  matrix  is  preserved  under  this  trans 
formation.  The  following  quantities  are  also  invariant  under  change  of  base: 

a.  det ([  S (x1  ,y ' )] )  *  det([T]^)  det  S(x,y)  det  T  ■  det  S(x,y) 


b.  Tr(P)  *  total  power  returned  to  a  pair  of  orthogonally  polarised 
antennas  *  Span([  S] )  ■  invariant. 

These  equations  are  quadratics  in  p,  where  p  relates  (x'y1)  to  (x,y).  For 
example,  if  [S(h,v)]  is  known  then  [S(left  circular,  right  circular)]  may  be 
predicted  by  setting  p  -  J.  Then 


[S(£,r)J 


shh"svv 
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SHH+SW 


SHH+SW 
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Substituting  the  values  for  the  trihedral  gives 


SU.r) 


n 


and  for  the  dihedral 


SU 


-  •  t: ;] 


This  confirms  the  fact  that  for  each  reflection  off  a  metallic  surface,  the  sense 
of  polarisation  is  changed. 


The  congruent  nature  of  the  change  of  base  transformation  has  several 
important  implications  for  interpretation  of  polarimetric  data  and  in  particular 
the  following  theorems  are  important: 

1.  Theorem  1 :  Any  square  matrix  [  S]  can  be  transformed  to  diagonal  form 
by  unitary  matrices  M  and  Q  such  that 

MSQ  ■  D  ■  o.  6 . _  6.-  Dirac  delta  function 

x  i  lJ  xJ 


Proof 


,  . *T  *T  *T  *T  * 

[MSQ]  -  Q  S  aM  -  a. 


Now  premultiply 


*T  *T  *T  *T  *T  * 

Q  S  M  MSQ  -  Q  S  SQ  -  a. ct. 6±J 


Since  Q  is  unitary 


Thus  Q  is  the  unitary  matrix  which  diagonalises  the  Hermitian  product 
S  S  by  a  similarity  transformation. 
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By  post-mule if  lying  we  have 
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Thus  M  is  the  unitary  matrix  which  diagonalises  the  Hermitian  product 

*T 

SS  by  a  similarity  transformation. 

T 

If  S  is  assumed  to  be  symmetric  then  M  *  Q 
Proof 


Hence  we  arrive  at  the  very  important  result  that  under  the  change  of 
base  transformation  the  coherent  scattering  matrix  can  always  be 
diagonalised 


This  is  a  result  of  fundamental  importance  since  it  implies  the 
existence  of  polarisations  which  remain  unchanged  when  incident  upon 
the  target.  These  so-called  eigenpolarisations  form  the  basis  for  the 
processing  techniques  to  be  described  later. 

2.  Theorem  2:  The  power  matrix  P  is  diagonalised  by  a  similarity  trans¬ 
formation  using  the  same  unitary  matrix  Q  which  diagonalised  the 
scattering  matrix  by  a  congruent  transformation  and  the  eigenvalues 
of  P  are  real  and  equal  to  the  squared  magnitude  of  the  eigenvalues 
of  [S] 

T* 

Proof  P  ■  S  S  and  so  by  1.  the  theorem  is  proved. 
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MEASUREMENT  OF  THE  SCATTERING  MATRIX 

There  have  been  very  many  different  techniques  proposed  for  measuring  the 

(21) 

target  scattering  matrix  and  a  good  review  of  these  is  provided  by  Huynen 
In  this  section  only  the  basic  requirements  will  be  outlined  and  some  of  the 
limitations  and  problems  pointed  out. 

The  scattering  matrix  with  absolute  phase  may  be  measured  using  a  coherent 
radar  with  a  dual  channel  receiver  capable  of  simultaneously  receiving  two 
orthogonal  polarisations  and  providing  I  and  Q  outputs  for  each  channel 
(Figure  15).  The  transmitter  must  then  transmit  two  orthogonal  polarisations 
such  that  the  four  complex  elements  of  [  S]  may  be  measured.  Ideally,  the  whole 
matrix  should  be  measured  at  the  same  time  but  this  requires  some  form  of  coding 
for  the  two  orthogonal  transmitted  polarisations  so  that  the  matrix  elements  can 
be  separated  on  receive.  The  easiest  way  of  doing  this  is  to  time  multiplex  the 
transmit  waveform  such  that  in  a  pulsed  radar  the  transmit  polarisation  is 
changed  on  a  pulse  to  pulse  basis.  Thus  it  takes  2  PRI  to  measure  the  matrix 
but  isolation  of  the  matrix  elements  is  guaranteed.  Note  that  if  a  priori 
knowledge  is  available  about  the  targets  to  be  measured  then  it  may  be  possible 
to  reduce  these  requirements  but  in  general  this  will  not  be  so. 

The  two  main  components  needed  for  the  SM  radar  above  those  of  a  conventional 
fixed  polarisation  device  are  a  polariser  in  the  transmit  channel  and  an  ortho¬ 
mode  coupler  (OMT)  in  the  receiver.  The  latter  device  is  a  standard  microwave 
component  and  allows  the  separation  of  an  incident  elliptically  polarised  wave 
into  two  orthogonal  components.  The  polariser  may  take  many  forms  such  as  a  mech¬ 
anical  switch  between  vertically  and  horizontally  polarised  feed  channels  or  a 
single  channel  Faraday  rotation  device  with  a  ferrite  phase  shifter  for 
generating  elliptical  polarisations.  Whatever  the  configuration,  the  transmitter 
must  be  able  to  switch  alternate  polarisations  at  the  PRF  rate.  This  brings  into 
play  an  important  problem  in  the  measurement  of  scattering  matrix.  If  the  target 
moves  during  the  2  PRI  time-period  necessary  to  measure  the  matrix  then  the 
measurement  will  be  inaccurate  and  if  movement  is  too  severe  and  the  time  period 
too  long  then  phase  determination  may  be  impossible.  Thus  coherent  measurement 
of  the  matrix  is  reserved  usually  for  carefully  staged  range  measurements  or  for 
high  PRF  radars  where  the  target  motion  is  known  to  be  within  a  certain  band¬ 
width  less  than  the  PRF. 
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The  problem  may  be  overcome  by  measuring  the  relative  scattering  matrix, 
since  then,  absolute  movement  of  the  target  over  the  period  of  matrix  deter¬ 
mination  is  less  critical.  A  schematic  diagram  of  a  pulsed  radar  suitable  for 
measurement  of  the  relative  SM  is  shown  in  Figure  14  together  with  a  typical 
measurement  switching  programme.  The  essential  difference  between  this  method 
and  the  previous  one  is  in  the  measurement  of  phase  angle.  This  latter  technique 
does  not  require  a  coherent  radar  since  the  PSD  may  take  one  of  the  received 
channels  as  its  phase  reference  and  so  provide  a  measure  of  the  time  phase 
difference  between  the  two  orthogonal  polarisation  components.  The  PSD  must  be 
able  to  unambiguously  determine  angles  in  the  range 

-180°  <  $pSD  <  180° 

One  interesting  possibility  is  to  mix  bases,  ie  transmit  on  base  (x,y)  but  receive 
base  (x',yf).  If  the  relationship  between  these  bases  is  known  as 

£(x,y)  -  [T]  E(x\y’) 

then  the  measured  matrix  [  Q]  may  be  converted  into  the  same  base  by 

E(x',y')  -  [QJ  E(x,y) 

-  [  QJ  f  T]  E(x\y’) 

This  technique  may  be  of  use  when  there  are  limited  dynamic  range  problems  in 
the  radar  receiver. 

To  date,  measurement  has  been  considered  using  ideal  antennae  and  a  noise 
free  system.  In  order  to  consider  a  more  realistic  system  it  is  necessary  to  be 
able  to  quantify  errors  introduced  by  using  real  antennas  which  have  finite 
losses  and  cross  polar  isolation.  In  order  to  do  this  one  may  consider  a  trans¬ 
mission  matrix  for  the  antenna.  In  other  words  if  polarisation  vector  P  is 
requested  for  transmission,  what  will  actually  be  transmitted  is 

-  [  E]  P 


For  an  ideal  system  [  E]  will  be  the  identity  matrix  but  for  real  systems  will 
represent  a  distortion  of  ellipticity  At  and  a  rotation  of  the  plane  of  polarisation 
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A9.  If  losses  within  the  antenna  are  assumed  small  then  [  E]  will  be  of  the  f 
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where  en,  e22,  e21  are  complex  and  Re(en)  <  0,  Re.(c22)  <  0,  the  equality 
holding  for  zero  losses,  are  functions  of  position  within  the  antenna  beam 
and  as  one  deviates  further  frcm  boresight  these  errors  tend  to  increase. 

On  receive  the  matrix  [  E]  operates  on  the  desired  polarisation  Q,  so  the 
net  effect  on  measurement  of  [  S]  is  the  congruent  transformation 


[  S]  -  [  E]  S]  [  E] 


[ E]  -  unitary  for  lossless  antenna 


In  reality  losses  may  arise  from  heat  loss,  reflections  from  support  structures 
and  radomes.  In  this  case  [  E]  is  non— unitary  and  the  determinant  and  Trace  of 
[  S]  are  no  longer  invariant  under  the  above  transformation.  Using  the  above  form 
for  [ E]  and  putting 


[S] 


S11  S12 


LS21  ^22' 


yields 
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Notice  that  if  S^*  are  small  relative  to  the  copolar  terms  then  because 
r'12  e 2}  are  small,  the  copolar  elements  are  less  susceptible  to  error.  The  cross- 
polar  terms  however  are  more  inaccurate. 


In  order  to  try  and  calibrate  these  errors  in  a  system  then  at  least  two 
calibration  targets  are  required  with  known  scattering  matrices.  For  example 
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hence  can  get  r.^*  £2l’  relate  to  the  cross-polar  isolation  of  the  antenna. 

If 
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and  if  the  absolute  RCS  of  calibrator  is  known  the  e^,  may  be  determined. 

Two  targets  which  correspond  to  these  calibration  matrices  are  a  trihedral  and 
dihedral  with  its  seam  at  45°  to  the  radar  vertical.  If  these  calibrations  are 
made  and  the  scattering  matrix  measured  within  the  target's  scattering  centre 
decorrelation  time  reliable  SM  data  may  be  collected  and  made  available  for 
processing.  Note  that  if  the  Stokes  parameters  are  used  then  they  give  a  direct 
measure  of  the  polarisation  purity  of  the  received  wave  and  its  susceptibility 
to  noise. 


TARGET  PARAMETERS 


By  measurement  of  the  five  parameters  of  the  relative  phase  SM  in  some  base 
(x,y),  the  scattering  for  any  incident  polarisation,  given  by  the  loci  of  points 
on  the  surface  of  the  PoincarS  sphere,  may  be  predicted  by  the  transformation 
equations  derived  earlier.  The  question  arises  as  to  how  best  process  this  data 
in  a  manner  which  will  yield  information  relating  to  target  geometry  and  symmetry 
independent  of  which  base  is  chosen  for  measurement.  It  will  be  shown  in  this 
section  how  such  a  set  of  five  target  parameters  may  be  derived  and  how  they  may 
be  used  to  improve  radar  performance. 
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The  target  descriptors  are  based  on  knowledge  of  the  so-called  target 
characteristic  or  null  polarisations.  These  are  transmit/receive  polarisation 
pairs  which  result  in  zero  target  backscatter.  There  are  in  general  four  such 
polarisations  and  every  target,  however  complex,  has  such  a  set.  The  proof  of 
the  existence  of  these  polarisations  is  based  on  theorem  1  from  the  previous 
section,  namely  that  there  always  exist  two  polarisation  vectors  which  diagonalise 
the  symmetric  scattering  matrix. 

The  two  polarisations  corresponding  to  the  diagonalisation  of  [  S]  are  termed 
the  eigenpolarisations  or  cross-polar  nulls  (XPOL)  in  that  when  they  are  trans¬ 
mitted  they  remain  unaltered  on  reflection,  so  that  by  receiving  in  the  orthogonal 
channel  a  null  response  would  be  obtained.  Care  must  be  taken  when  calculating 
these  polarisations  since  as  usual  the  co-ordinate  frames  must  be  related  so 
that  the  Poincar£  sphere  representation  may  be  used.  When  allowance  is  made  for 
the  conjugate  nature  of  the  backscatter  co-ordinates  when  compared  to  the  trans¬ 
mitter  then  the  eigenpolarisations  are  solutions  of 


IS)  S* 


A  -  complex  eigenvalues 


It  is  now  apparent  why  these  are  called  eigenpolarisations  since  this  equation 
is  similar  to  the  classical  eigenvalue  problem 

[S]  £  -  A£ 

The  conjugate  sign  in  the  polarisation  problem  forces  the  corresponding  eigenvalue 

to  be  phase  determined  whereas  in  the  more  commonly  met  form  of  eigenvalue  problem, 

o 

if  X  is  a  solution  then  so  is  Ae  0  <  6  <  360  . 

From  theorem  2  the  eigenvectors  of  [  S]  are  the  same  as  those  of  S*S  and  the 
latter  has  real  eignvalues  given  by  the  square  moduli  of  the  eigenvalues  of  [  S] . 
When  the  eigenvalues  of  [  S]  are  distinct  then  it  may  be  shown  that  the  two 
corresponding  eigenvectors  are  always  orthogonal. 


Proof  let  Aj,  A^  be  the  two  eigenvalues  and  A^  >  A^ 
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Since  [ S]  is  symmetric 

l  S] £x  *  £2  -  £x  •  [S}£2 

l*J  l£j*  •  £2I  “  lx2(  1^1  *  h* 


hence  if  X.,  ¥>  X 2 
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This  is  the  condition  for  othogonality  of  polarisation  vectors. 

On  the  Poincar£  sphere  £^«£2  lie  diametrically  opposite  and  so  knowledge  of 
one  immediately  determines  the  other.  Thus  two  of  our  desired  target  parameters 
are  the  latitude  and  longitude  or  ellipticity  and  inclination  angle  of  one  of  the 
eigenpolarisations.  The  one  chosen  is  the  eigenvector  with  the  largest  corres¬ 
ponding  eigenvalue  since  the  maximum  RCS  of  the  target  is  given  by  the  square 
modulus  of  the  maximum  eigenvalue. 


Proof 
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Hence  if  |x^|  >  |X2|  then  maximum  power  is  obtained  if  |b| 
polarisation  E  ■  E,  and  P  *  i X , I2.  In  the  above 
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This  maximum  polarisation  is  the  third  target  parametar  and  is  obviously  of 
great  importance  to  radar  systems  since  it  is  the  maximum  RCS  that  may  be 
obtained  for  a  target  at  a  fixed  aspect  and  frequency. 

Another  method  of  calculating  the  eigenpolarisations  is  to  consider  the 
general  transformation  equations,  quadratic  in  p.  The  eigenpolarisations  are 
seen  to  be  solutions  to 

SXY  “  0  pSYY  ~  P*SXX  +  SXY(1  “  pp  ) 

it 

This  equation  may  be  solved  by  noting  that  S  S  has  the  same  eigenvectors  as  [  S) 
and  is  diagonalited  by  a  similarity  transformation 

T*1  S*S  T 


If  this  matrix  product  is  expanded  then  for  the  off-diagonal  terms  to  be  zero 
the  equation 

p2b  +  p(a-c)  -  b*  -  0 


where 


3  "  SXXSXX  T  SXYSXY 


b  “  SXYSXX  +  SYYSXY* 


c  -  S^S^*  +  SxySxY* 


so 
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These  equations  give  two  values  ox  p  corresponding  to  the  two  eigenpolarisations. 
By  writing  p  *  tana  e^  then  the  latitude  and  longitude  of  the  eigenpolarisations 
may  be  found  from 


0^  ■  j  arctan(tan2a  coefi) 


■  J  arcsin(sin2a  sini) 
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As  a  simple  example  of  the  eigenpolarisations  of  a  target  consider  the  back- 
scatter  from  a  wire  grid  target  as  shown  in  Figure  16.  If  the  grid  is  aligned 
at  0°  to  the  incident  vertical  then  the  polarisation  of  the  backscattered  wave 
will  be  vertical,  independent  of  the  incident  polarisation.  Thus  vertical  is 
one  eigenpolarisation  while  horizontal  is  the  other.  Further,  vertical  is 
obviously  the  maximum  polarisation  since  for  horizontally  polarised  waves  there 
is  zero  backscatter.  If  the  grid  is  now  rotated  about  the  radar  line  of  sight 
then  these  eigenpolarisations  remain  linear  but  rotate  in  the  equatorial  plane 
of  the  Poincard  sphere. 

The  remaining  two  characteristic  polarisations  are  the  copolar  nulls  (COPOL) 
or  incident  polarisations  which  are  transformed  into  their  orthogonal  state  on 
reflection.  These  may  be  determined  from  the  transformation  equations  by 


sxx  -  0 


p  Sy„  +  2S^p  + 


"SXY  1  *^XY  "  SXXSYY 


It  is  these  polarisations  which  are  exploited  in  rain  clutter  suppression  radars 
since  spherical  raindrops  have  left  and  right  circular  polarisation  as  their 
copolar  nulls.  Notice  that, although  in  this  case  the  COPOL  nulls  are  orthogonal, 
in  general  they  will  not  be  so. 


Of  particular  importance  is  the  relation  of  the  COPOL  nulls  to  the  eigen¬ 
polarisations  and  in  order  to  develop  this  relationship  consider  transformation 
of  the  scattering  matrix  when  it  is  expressed  in  (X^^)  where  X^  and  X2  are  the 
eigenpolarisations  , , 
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For  the  transformed  matrix  to  be  a  COPOL  null  then  o  must  satisfy 
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From  this  we  may  write  [  D]  in  the  form 


J(*+v) 


.  2  J(<p-v) 

m  tan  ye 


where  m  ■  A, 


Y  -  90°  -  o  0  <  y  <  45° 


The  two  solutions  satisfying  these  conditions  are 


v  -  6  -  90°  -90°  <  v  <  90° 


.  ,  -  /  1  „  J 6 

P3  +J  “  tana  e 
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In  terms  of  the  Deschamps  co-ordinates  on  the  Poincarl  sphere  these  points  are 
given  by  (2yj-v)  and  (-2Yrv)  taken  from  the  minimum  eigenpolarisation  X2-  In 

A 

order  to  express  the  transformations  from  the  maximum  polarisation  X^,  then 
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2y  ■*■  180°  -  2y  and  v  *►  -v.  Hence  the  points  become  (2a, v)  end  (-2a, v),  where 
positive  v  is  taken  in  an  anticlockwise  direction  from  the  equatorial  plane  and 
Y,v  are  angles  from  ID]  . 


These  results  are  very  important  since  they  show  an  interesting  geometrical 
relationship  between  the  two  sets  of  null  polarisations.  As  has  already  been 
noted  the  eigenpolarisations  lie  diametrically  opposite  and  now  we  have  the 
additional  information  that  the  C0P0L  nulls  have  co-ordinates  (±  2a,  ~v)  from 
the  maximum  eigenpolarisation.  Thus  all  four  points  lie  on  one  great  circle 
which  defines  a  plane  in  polarisation  space  (Figure  17)  and  the  diameter  joining 
the  eigenpolarisations  bisects  the  planar  angle  between  the  COPOL  nulls.  The 
resulting  structure  is  known  as  the  polarisation  fork  and  is  a  powerful  graphical 
aid  to  radar  target  classification  studies.  The  fork  prongs  are  drawn  from  the 
origin  to  the  COPOL  nulls  and  the  length  of  the  'handle'  is  the  radius  of  the 
Poincard  sphere  or  the  span  ([  S] ),  which  is  a  transformation  invariant  equal  to 
the  power  in  the  polarised  component  of  the  wave.  Note  that  some  authors  take 
the  target  maximum  RCS  as  radius  of  the  sphere  but  in  order  that  the  results 
be  consistent  with  the  Stokes  formalism  it  will  be  taken  in  these  lectures  as 

r  “  812  +  g22  +  s32  "  Sp(t  SJ  > 

This  now  enables  us  to  classify  different  radar  targets  according  to  a  five 
element  target  vector  given  by 

V  -  (20J,  2t ^ ,  v,  y ,  m) 

The  first  three  angles  define  the  plane  of  the  fork  in  three  dimensions,  the 
fourth  angle  is  y  being  the  fork  separation  and  relating  the  difference  in 
magnitude  of  target  eigenvalues.  The  final  parameter  m  is  the  target  maximum 
RCS.  In  the  next  section  these  ideas  will  be  developed  further  and  examples 
given  but  first  consider  again  the  diagonalised  form  of  [  S] 
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Using  tha  Pauli  matrix  algabra  developed  earlier  we  may  write 
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This  form  is  used  when  developing  relationships  between  the  coherent  theory  and 
Stokes  matrix.  This  is  the  next  logical  step  in  the  theory  but  is  beyond  the 
scope  of  these  lectures.  For  those  interested,  further  details  may  be  obtained 
in  references  5  and  15.  (See  also  appendix  1). 
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CLASSIFICATION  OF  RADAR  TARGETS 


In  the  previous  section  the  target  vector  V  was  developed  which  has  five  j 

parameters  end  presents  a  full  polarimetric  descriptor  of  the  target.  In  this 
chapter  some  time  will  be  spent  investigating  the  form  of  V  for  some  simple 
radar  ocatterers  and  outlining  a  basic  scheme  for  classifying  radar  targets 
according  to  their  polarimetric  properties.  Note  that  we  are  still  considering  1 

fixed  aspect  and  frequency  data  and  coherent  scattering  in  the  far  field  of  the 

target.  The  problems  of  scintillating  targets  and  change  of  aspect  and  frequency  \ 
will  be  briefly  considered  later.  1 

i 

>i 

» 

Let  us  define  the  polarisation  spectrum  of  a  target  with  scattering  matrix  < 

[  S]  by  the  copolar  and  crosspolar  RCS  both  given  as  a  function  of  polarisation.  ■! 

j 

Shown  in  Figures  18  and  19  is  the  spectrum  of  a  simple  flat  plate  or  trihedral  y 

target.  On  these  plots  6  is  mapped  from  0  ■+  180°  versus  t  the  ellipticity  J 

(~45°  <  t  <  45°).  Notice  that  the  copolar  plot  shows  a  maximum  RCS  for  all  J 

linear  polarisations  and  a  null  for  left,  and  right  circular.  The  crosspolar  «! 


plot  shows  a  maximum  for  circular  and  a  corresponding  null  for  all  linear. 

These  properties  may  be  summarised  by  writing  the  target  vector 

V  ■  (0°,  0°,  0°,  45°,  m)  m  ■  A  -  area  of  plate 

X  -  illuminating  wavelength 

When  mapped  onto  the  Poincar£  sphere  the  polarisation  fork  for  this  target  is  as 
shown  in  Figure  20.  The  fork  has  maximum  angle  separation  (y)  which  is  indicative 
of  the  copolar  nulls  being  orthogonal.  Notice  that  because  the  eigenpolar isations 
include  the  whole  equator  this  target  is  rotation  invariant  ie  there  is  no  change 
in  the  fork  if  the  target  is  rotated  about  the  radar  line  of  sight. 

Compare  these  results  with  those  for  the  dihedral  reflector  as  shown  in 
Figures  20-23.  The  dihedral  is  assumed  to  have  its  seam  parallel  to  the  radar 
vertical  and  thus  has  two  nulls  at  ±  45°  linear  and  maxima  for  all  polarisations 
with  0  ■  0°  or  90°  and  -45°  <  t  <  45°.  The  target  vector  is 

V  -  (0°,  0°,  90°,  45°,  m) 

The  polarisation  fork  is  similar  to  that  of  the  trihedral  except  for  a  90° 
rotation  about  the  eigenpolarisation  axis.  This  is  the  skip  angle  and  relates 
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Co  the  number  of  bouncas  tha  reflected  aignal  undergoes.  In  general  however  the 
skip  angle  is  the  phase  angle  between  eigenvalues  and  so  does  not  lend  itself  to 
such  simple  physical  interpretation.  Note  that  because  of  the  form  of  polarisa¬ 
tion  fork  the  dihedral  is  rotation  dependent  and  for  a  general  rotation  angle 
the  target  vector  becomes 

V  -  (2tp,  0,  90°,  45°,  m) 

ie  the  fork  is  rotated  in  the  equatorial  plane  by  2t />° . 

The  wire  grid  target  considered  earlier  has  the  polarisation  fork  shown  in 
Figure  24.  In  this  case 


V  -  (2<K  0°,  0°,  0°,  m) 


I 


f 

1 


Notice  that  the  copolar  nulls  are  now  coincident  and  the  fork  has  collapsed  into 
a  straight  line.  For  the  grid  at  an  arbitrary  angle  \p  it  matrix  is 


t  S]  - 


cos  y  Jsin2t/i' 


2 

,Jsin2^  sin  4>. 


and  as  such  the  fork  rotates  by  2\p  in  the  equatorial  plane. 


These  are  all  examples  of  an  important  class  of  symmetric  targets.  The 
general  form  of  [  S]  in  (x,y)  for  a  symmetric  target  is 
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re2iv 
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tan  y 


-2iv 


where  (x,y)  is  a  linear  base 


All  targets  having  an  axis  of  roll  symmetry  are  symmetric  at  any  aspect  angle, 
eg  cones,  cylinders,  ellipsoids  and  combinations  of  these.  Other  targets  like 
corner  reflectors  may  have  a  plane  of  symmetry  through  a  line  of  sight  direction 
and  thus  be  symmetrical  targets.  All  these  scatterers  have  a  target  vector 

V  *  (6,  0°,  v,  y,  m)  ie  t  *  0° 
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This  may  ba  illustrated  by  virtue  of  Figure  25  which  shows  a  general  symmetric 
radar  target  and  a  proposed  maximum  polarisation  ?inax*  Because  of  the  symmetry 

A 

of  the  object  there  must  then  be  another  maximum  polarisation  P  obtained 

max  z 

from  P  by  reflection  in  the  plane  of  symmetry.  However  theory  states  that  the 

A 

eigenpolarisations  are  orthogonal  and  so  Pmax  must  be  either  parallel  or  perpen¬ 
dicular  to  the  symmetry  axis.  For  either  choice  Tm  "  0  and  the  eigenpolarisations 
are  linear. 

Symmetric  targets  are  characterised  by  their  eigenpolarisations  always  lying 
in  the  equatorial  plane  and  are  of  special  interest  in  radar  scattering  studies 
since  many  of  the  most  common  radar  targets  show  these  features. 

The  class  of  non-symmetric  targets  are  classified  by  *=  0.  The  maximum 

value  that  t  may  achieve  13  ±  45°  and  it  does  this  for  a  helix  target  which 
m 

has  the  matrix 

J  “  Tl  -J' 

or 

^left  screw  *“  ^  bright  screw 

The  corresponding  forks  are  shown  in  Figures  26  and  27  and  for  these  targets 
V  -  (0°,  ±  90°,  0°,  m) 

Because  these  represent  an  extreme  in  target  asymmetry  the  angle  Tm  is  termed 
the  helicity  angle. 

A  class  of  nonsymmetric  targets  which  has  special  significance  in  the 
decomposition  theorems  of  Huynen  are  n-targets  (nonsymmetric  noise)  defined  by 

t  *  ±  45° 

m 

Thus  [  S]  has  the  general  form 
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Note  that  the  dihedral  target  may  be  considered  as  a  special  case  of  an  n-target 
since 


[S] 


cos2a 

sin2a 


sin2a 

-cos2a 


So  far  only  single,  isolated  targets  have  been  considered  but  in  reality  radar 
targets  will  be  composed  of  many  scattering  centres  each  with  different  polari- 
metric  properties  and  target/radar  motion  will  result  in  scintillation  of  the 
target  echo.  In  the  fixed  polarisation  scalar  theory  these  fluctuations  have 
received  much  attention  but  little  is  known  of  the  dynamic  properties  of  the 
target  vector  V  and  their  relation  to  changes  of  radar  frequency  and  target 
aspect.  Certain  interesting  trends  may  be  pointed  out  however. 


If  the  target  has  a  few  dominant  scatterers  with  different  polarimetric 
properties  then  the  loci  of  null  polarisations  for  charge  of  aspect  or  frequency 
might  be  expected  to  cover  large  areas  of  the  Poincarl  sphere  but  if  on  the  other 
hand  it  consists  of  reflectors  with  similar  polarimetric  properties  or  many 
scatterers  distributed  within  a  narrow  dynamic  range,  then  clustering  of  the 
null  polarisations  might  be  expected.  This  has  been  verified  for  rain  clutter^  ^ 
and  provides  the  basis  for  polarisation  filtering  whereby  controlled  suppression 
may  be  had  for  limited  areas  of  polarisation  space.  The  frequency  dependence  of 
V  is  an  important  area  of  study  and  has  been  covered  by  Boerner  .  This  work 

promises  to  combine  polarimetric  information  with  spatial  range  information 
available  from  broadband  systems  and  as  such  provide  a  powerful  theory  for 
electromagnetic  inverse  scattering. 


The  incorporation  of  full  polarimetric  information  with  radar  systems  has 
many  widespread  applications.  The  efficient  suppression  of  rain  clutter  is  but 
one  example  and  in  general  any  target  whose  nulls  show  clustering  in  polarisa¬ 
tion  space  may  be  enhanced  or  suppressed  using  these  techniques.  The  important 
point  to  note  is  that  SM  methods  allow  the  choice  of  effective  transmit/receive 
polarisation  to  be  made  in  the  signal  processing  and  real  time  adaptation  of 
these  is  then  limited  only  by  the  speed  of  processing  available. 


The  techniques  have  application  in  improving  signal  to  clutter  ratio,  ECCM, 
multipath  reduction  and  radar  target  identification  studies.  In  particular  the 


Sys  are  the  coefficients  of  expansion  for  the  transformed  coherency  matrix. 
Writing  this  in  compact  form  and  using  the  cyclic  property  of  the  Trace  operator 
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where 
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The  coherency  matrix  may  be  written 
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We  have  now  found  a  matrix  M  whose  element  Myv  relates  the  vth  Stokes  parameter 
of  input  wave  to  the  yth  Stokes  parameter  of  reflected  wave  where 
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Thi(°  relationship  gives  a  method  of  determining  the  Stokes  reflection  matrix  M 
from  any  given  target  matrix  [  T] 
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which  after  taking  the  trace  gives 
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This  expression  for  l  M]  is  in  the  natural  ordering  and  in  order  to  convert  to 
traditional 
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Thus,  using  the  same  notation  as  above 
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Before  working  out  an  explicit  form  for  [ Mj  care  must  again  be  taken  to  ensure 
that  the  incident  and  reflected  waves  are  referenced  to  the  same  co-ordinate 
system,  f M]  has  already  been  corrected  for  x,y  plane  co-ordinates  through  [  T] 
but  in  order  that  the  conjugate  nature  of  reflected  wave  co-ordinates  be  taken 
into  account. 


0  0 

0  0 

1  0 

0  -1 


[QUe 


"  l  QH  M]  & 


If  we  then  define 


[Mj  -  1(01  HI 
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Chen  we  may  use  ( M]  '  Co  relaCe  an  incident  wave  P  to  a  reflected  wave  where 
both  are  defined  in  cerms  of  the  transmitter  co-ordinates  and  write  the  expression 
for  received  power  P  as  a  dot  product  of  Stokes  vectors 


P  -  [  M3  '  £(P)  g(Q) 


£(P)  -  Stokes  vector  for  P 
£(Q)  -  Stokes  vector  for  Q 


When  multiplied  out 
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then  [  M] '  becomes 
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Notice  that  for  symmetric  [  Tj ,  [  M]  '  is  also  symmetric  with  nine  independent 
elements.  For  the  general  bistatic  case  however  there  will  be  sixteen  elements 
although  symmetry  in  the  scatterer  may  reduce  this  number. 
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For  partially  polarised  waves,  averages  are  taken  of  the  individual  Stokes 
parameters  and  so  one  may  define  an  average  Stokes  reflection  matrix  by 

IM]  '  ■  <[  Ml  >  where  m.  ,  ■  <m.  ,> 

1  ^ov  1  iJav  lJ 

Huynen's^J  decomposition  is  based  on  the  fact  that  this  averaged  matrix  may  be 
decomposed  into  a  component  with  corresponding  coherent  matrix  [  T]  satisfying 
the  condition  for  physical  realisability  of  Stokes  vectors,  plus  a  remainder 
matrix  which  may  be  decomposed  itself  into  the  sum  of  two  n-targets. 
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FIG.  2  ORTHOGONAL  POLARISATIONS 
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FIG.  7  ANTENNA  CO-ORDINATE  SYSTEMS 
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